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Abstract. The conditions of thermodynamic equiHbrium, local ther- 
modynamic equilibrium, and statistical equilibrium are discussed in 
C/j ■ detail. The equations of statistical equilibrium and the supplementary 

(-H ' equations are shown together with the expressions for radiative and 

O ,. coUisional rates with the emphasize on the solution for trace elements. 

o 

H 

^ 1 Introduction: Stellar atmosphere problem 

Solution of the stellar atmosphere problem is usually considered as th e deter 



ruination of spatial dep endence of basic macroscopic quantities (see also iKubat , 



I> ' bOldlthese proceedingsf ) also be understood (following iHubenvl . 1 1 9 7d ) as aiming at 



OO finding three basic microscopic distributions, namely, the momentum distribution 
^jj (distribution of velocities of all particles), distribution of particle internal degrees 
p.^ , of freedom (populations of atomic excitation stages), and distribution of internal 

degrees of freedom of the electromagnetic field (radiation field for all frequencies 

and directions). 



o 
o 



^ 



2 From thermodynamic equilibrium to NLTE 

2.1 Thermodynamic equilibrium 



^ ,. To main tain thermody na mic equilibrium , several necessary conditions have to be 
■ ■ ' met (see iHubenvll 1 9 76l also lHubenvlll997 ). First, the characteristic time for macro- 



met (see 

scopic changes of the medium has to be much larger than the relaxation time 
(^relaxation < ^macroscopic ciianges)- Sccoud, macroscopic chaugcs of the medium 
must be small on the scale of a particle mean free path (^macroscopic ciiangcs ^ 
'free patii)- Then, the relaxation time has to be much less than the time inter- 
val between inelastic colhsions (treiaxation < ^inelastic collisions)- And, finally, if the 
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latter condition is violated, then colliding particles need to have equilibrium dis- 
tributions. 

If the medium is in thermodynamic equilibrium, then all three distributions 
mentioned in Introduction (momentum, internal, and radiation) have their equi- 
librium form. Velocities Va of particles a follow the Maxwellian distribution, 

/(«a) d«a = ^e-«/<) dv,, (2.1) 

where Ua = y'2kJ]Jrri^ is their most probable speed, Ta is the corresponding 
particle temperature, rua, is its mass, and k is the Boltzmann constant. In ther- 
modynamic equilibrium all temperatures Ta are the same and in the following we 
denote them as T. 

Individual at omic energy le vels i are populated according to the Boltzmann 
distribution fsee lMihalal Il978l Eq. 5-4), 



± -_ 


= ^e-^ 





9o 



(2.2) 

where Ui is the population of the i-th level, no is the population of the ground 
level, asterisks (*) denote equilibrium values, gi is the statistical weight of the level 
i, and Xi is its exc itation energy. The ionization equilibrium is described by the 



Saha equation (see iMihalasl . Il978l Eq. 5-15) 



n; ^ u,{T) ( 



7V;+i 2Uj+i{T) \2TTmekT 



e^, (2.3) 



where N* is the total population (number density) of the ion j, Uj{T) is its 
partition function, n^ is the electron number density, xij is the ground state 
ionization energy, and h is the Planck constant. 

Finally, the radiation field energy is distributed in frequency v according to 
the Planck law 

BAT)^'-^^. (2.4) 

C e^r — 1 

2.2 Local thermodynamic equilibrium (LTE) 

The fact that we can observe stellar light and the shape of the stellar spectrum 
clearly show that the stellar radiation is far from that of a blackbody. Therefore, 
the assumption of radiation field energy equilibrium is not valid and we have to 
include detailed radiative transfer. However, simple relaxing of the equilibrium 
condition for radiation does not help, since if the first two distributions are in 
equilibrium in a homogeneous medium, then also the radiation field is forced to 
have the equilibrium value. 

To allow the light to escape from the star, we have to allow spatial dependence 
of temperature and density. We assume local equilibrium distributions of particle 
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velocities (12.11) and level populations (|2.2l and l2.3|) . which means that the equilib- 
rium condition ignores possible gradients of T{f) and N{r). Then, the equilibrium 
conditions change from point to point. This approximation is usually called local 
thermodynamic equilibrium (LTE). 

To determine the specific radiation field intensity (Ifn,) distribution, we solve 
the radiative transfer equation. For the simple case of a plane-parallel horizontally 
homogeneous atmosphere (ID case with the only spatial coordinate z ), the transfer 
equation along the ray determined by /x = cos 9, where 9 (see Fig. 1 in lKubati . l20ld 



these proceedings! ) is the ray angle with the coordinate axis z, reads 



dz 



M-J7^ = ?7i/ - Xi^V' (2-5) 



where 7]^, and Xv ^'^^ the emissivity and opacity, respectively. The condition of 
local thermodynamic equilibrium implies that 77,^ == X'jB^{T{z)) at each depth z. 

2.3 Statistical equilibrium (NLTE) 

The non-equilibrium distribution of the radiation field obtained by solution of (J2.5I) 
affects in turn the atomic level population distribution. If this influence becomes 
strong, the Saha-Boltzmann ionization and excitation distribution is no longer in 
equilibrium. But even if this influence is weak, it is better to assume that the 
excitation distribution can generally be a non-equilibrium one. 

Doing this, we assume that the excitation and ionization states distribution 
does not follow from the Saha (12. 3p and Boltzmann (j2.2[) equations. For its de- 
termination we have to explicitly consider all individual level populating and de- 
populating processes. Such treatment is usually referred to as NLTE or non-LTE. 

Assuming NLTE we still assume an equilibrium (i.e. Maxwellian - Eq. 12. 1[) 
distribution of particle velocities, especially of electrons. The radiation field and 
level populations are allowed to have non-equilibrium distributions, which for the 
radiation field is determined by solution of the radiative transfer equation (|2.5I) 
and for the level populations by solution of statistical equilibrium equations (|4.3p 
described later in the Section ID 

2.3.1 Microscopic processes 

Due to the non-equilibrium nature of level population distribution it is necessary 
to take into account all processes infiuencing level populations on a microscopic 
level. The most important processes for hot star atmospheres are listed in Table [H 

2.4 Between LTE and NLTE 

2.4.1 Elastic collisions 

In both LTE and NLTE, a Maxwellian velocity distribution of particles (|2.ip is 
assumed however it is never explicitly stated. Elastic collisions between parti- 
cles tend to maintain their equilibrium velocity distribution. These collisions are 
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Table 1. List of most important microscopic processes affecting the equilibria in stellar 
atmospheres. In this table, v means photon, e stands for electron, X means a particular 
atom. Superscript ^ denotes the ion and * the excited state. E (in the rightmost column) 
means that the process tends to establish equilibrium of population numbers, M denotes 
processes that tend to establish equilibrium velocity distributions. 

collisions 

elastic (e-e, e-H, e-H+, e-He, H-H, H-He, ...) M 



inelastic with electrons 
excitation 
de-excitation 
ionization 
recombination 



e(w)+X^e(w' <?;)+X* E 

e(w)+X*^e(w' > w)+X E 

e(w) + X ^ e(w' < w) + e(i/') + X+ E 

e + e(w)+X+ ^e(u' 7^w)+X E 



inelastic with other particles less frequent => neglected E 


interaction with radiation 
excitation 




v + X^X* 




de-excitation 


spontaneous 
stimulated 




X* ^v + X 

v + X* ^Iv^X 




ionization 


autoionization 
Auger ionization 




1/ -H X ^ X+ + e 

J/ -K X ^ X** ^ X+ -h 6 

J/ -F X ^ X+* + e 




recombination 


spontaneous 

stimulated 

dielectronic recombination 


e -h X+ ^ !/ + X 

V ^c + X^ ^2v + X 

X+ + e ^ X** ^ 1/ + X 


E 
E 
E 


free-free transitions 




zy + c + Xf^ c + X E 




electron scattering- 


free (Compton 
bound (Raylcig 


Thomson) i^ + c — > i^ + e 
;h) v + X^v^X 





rarely explicitly considered both in LTE and in NLTE calculations. Usually, it is 
sufficient to assume that they are frequent enough to balance all non-equilibrium 
processes that tend to destroy the Maxwellian velocity distribution. For stellar 
atmospheres, t he equilibrium velocity distribution is a reasonable assumption (see 
Mihalasl . [l978h . 



2.4.2 Inelastic collisions 

Inelastic collisions play a key equilibrium role in the atomic level populations. If the 
velocity distribution is in equilibrium, then the inelastic collisions push the atomic 
level populations also to their equilibrium. On the other hand, inelastic collisions 
change the electron kinetic energy and so cause departures from the Maxwellian 
velocity distributions. Equilibrium velocity distribution is recovered only if there 
is enough elastic collisions. Luckily, in most astrophysical cases the time between 
two inelastic collisions is much larger than the time necessary for rebuilding the 
equilibrium velocity distribution - the relaxation time (the inequality ^relaxation <C 
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^inelastic collisions tiolds), SO we cau safsly adopt Maxwellian velocity distribution. 
However, there are exceptions, for example there is evidence that the electron 
velocity distr ibution differs from the Maxwellian in t he solar transition region 



ity 



( Shoubl . 119831 ) and in the solar corona (jCranmeii I1998I ). In such cases we would 



have to solve an additional kinetic equation, namely the equation for electrons. 

2.4.3 Level populations 

In the following, we assume a Maxwellian velocity distribution for all particles. 
The distribution of the radiation field is, on the other hand, not in equilibrium 
and must be determined by solution of the radiative transfer equation. Whether 
level populations will have their equilibrium value or not (in other words, if the 
atmosphere will be "LTE" or "NLTE" ) depends on the balance between processes 
connected with the equilibrium part - inelastic collisions with electrons having 
Maxwellian distribution - and processes connected with the non-equilibrium part 
- radiative processes, which connect the level populations with a generally non- 
equilibrium radiation field. In Table [T] the letter (E) denotes processes that aim 
at maintaining equilibrium of level populations. 

If the equilibrium processes dominate (e.g. in a coUisionally dominated at- 
mosphere), then level populations will have or will be close to the equilibrium 
Saha-Boltzmann values. In such case, all processes are in detailed balance, which 
means that the rate of each process is balanced by the rate of the reverse pro- 
cess. For a Maxwellian velocity distribution, the coUisional processes are always 
in detailed balance. Note that all radiative processes involving a free electron (re- 
combination, free-free transition) are in principle collisional. On the other hand, 
radiative processes can hardly be in detailed balance, they may be in detailed bal- 
ance only for the case when radiation has its equilibrium, i.e. Planck, distribution, 
which happens only in the deepest parts of stellar atmospheres. When the radia- 
tion intensity is different from Planckian, J^ 7^ B^, then LTE is not an acceptable 
approximation. 

If the non-equilibrium processes dominate, then a quite complicated game be- 
gins. The medium may reach states very far from equilibrium. 



3 Trace elements 

Solution of the NLTE problem for trace elements may be considered as a sim- 
plification of the standard model atmosphere problem, where we aim to deter- 
mine the tempera ture and density structures of the atmosphere (see lKubatl . I2OIOI 



these proceedinggi ). Here we assume that we have a given model atmosphere (ei- 



ther LTE or NLTE), i.e. given temperature {T{r)) and density (ne(r)) structures, 
given chemical composition, and if the given model atmosphere is a NLTE one, 
also occupation numbers of atomic levels for the most important atmosphere con- 
stituents. This implies that also the background opacities (which are calculated 
from the given model atmosphere) are known. Then we solve simultaneously the 
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radiative transfer and statistical equilibrium equations for a trace element (or trace 
elements) . 



There are several conditions, which have to be met to enable the element to be 
a trace element. The basic requirement is that the influence of the trace element 
on the atmospheric structure is really negligible. The trace element may, of course, 
have influence on emergent radiation in some parts of the spectrum. More exactly, 
there is a requirement that inclusion of the NLTE trace element does not change 
the atmospheric structure, which may also allow to include the case of a change 
from an LTE treatment of the element to the NLTE one. This may sometimes 
enable to conside r even hydrogen as a trace element (see also IPrzvbillal . uOlOl 



these proceedings! ) ■ but this is not generally recommended. Stellar atmospheric 



structure is usually sensitive to changes in the dominant ions, even if such changes 
are relatively small. Sometimes it is very difficult to predict the final effect caused 
by a change of the basic assumption from LTE to NLTE, so after any NLTE trace 
element calculation it is necessary to check if the assumptions we made at the 
beginning are still met. 



This also means that we have always to check if a trace element is really a trace 
element. For example, changes in ionization equilibrium alter also the free electron 
density, which may have an impact on the ionization structure of the initial model 
atmosphere. Background opacities have to be the same in the solution of the trace 
element as in the model atmosphere. 



There is one principal warning for using background LTE models, which is 
based on physical grounds. As a commonly used case we have an LTE model 
atmosphere, it means that we assume that in the whole atmosphere there are 
sufficient collisions with electrons to maintain equilibrium values of population 
numbers. If we then choose one of elements as a trace element and solve the 
NLTE problem, we implicitly say that the equilibrium may be violated for this 
particular element. But since the number of free electrons is the same as for all 
other elements (which populations have their LTE values thanks to the assumption 
of the LTE background model atmosphere), we simultaneously assume that for 
the same physical conditions there exists detailed balance for all transitions that 
belong to non-trace elements, but at the same time it is lacking for the trace 
element. Thus we are using two contradictory physical assumptions for the same 
region. Since LTE is a special case of a more general NLTE case, there may be 
regions in the atmosphere, where this very specific condition is fulfilled, but it 
hardly happens for the whole atmosphere. Using an LTE model atmosphere and 
assuming NLTE for a trace element is highly inconsistent and one has always to 
prefer NLTE model atmospheres as background models. 
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4 Equations of statistical equilibrium 

For the change of the occupation number n, of t he state i of each element in the 



stellar atmosphere we may write ([Mihalasl . I1978L Eq. 5-48) 



o NLtot 

-^ + V • in,v) = J2 i^iPi^ - ^^Pa) ' (4-1) 

1=1 

where NLtot is the total number of energy levels considered, Pu is the probability 
(per time unit) of a transition from the state (level) i to the state (level) I. This 
probability consist of two parts, namely of the probability that the transition is 
caused by absorption or emission of radiation (the radiative probability Ru), or 
that the transition is caused by a collision with a neighbouring particle, mostly 
with an electron (the coUisional probability Cu). Then, the probability Pu can be 
expressed as a sum of radiative and collisional probabilities, Pu = Ru + Cu. If 
these probabilities are multiplied with the starting level population, then we call 
them radiative and collisional rates, respectively. 

For the case of a stationary atmosphere, i.e. if the time variations of the 
medium are negligible, we may neglect the time derivative in the equation (|4.1|) 
and we obtain the time independent set of the statistical equilibrium equations, 

V■in^v)^J2 ^'^iPi^ - "^^Pa ) ■ (4.2) 

This set of equations has to be used in stationary moving atmospheres, which 
includes the common case of a stellar wind. However, due to simplicity, the ad- 
vective term is very often neglected. As a result we obtain an equation, which is 
valid in static atmospheres exactly, 

Y,{nlPl^-n^P^l)=0. (4.3) 

We say that the particular transition i <-> / is in detailed balance, if niPu — riiPu- 
If this condition is fulfilled for all j, I, then all transitions are in detailed balance, 
which is the case of LTE. Then, also the level populations have their equilibriiun 
values n* given by Eq. (|2.2p . 

4.1 Equilibrium level populations 

Equilibrium level populations are calculated using the assumption of thermody- 
namic equilibrium or local thermodynamic equilibri um. They are u sually denoted 
as n* and referred to as LTE populations. However. iMihalasI ( 19781 ) uses the quan- 



tity n* in a different meaning. The LTE populations are calculated wi th respect 
to th e ground level of the next higher ion using the equation (Eq.5-14 in iMihalasl . 
19781) 

3. 

'■^ 50J+1 2 \2TTmekT J 
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Here j is the index of the corresponding ion and noj+i is the actual population 
of the ground level of the next higher ion, which is not necessarily an LTE value. 
Equation (I4.4p is used as a definition of LTE populations in NLTE c alcula tions. 
Then, if we are using the departure coefficients b^ = rii/n* ( Menzei Il937 ). it is 
necessary to know which definition of LTE population is used. Note that for both 
definitions of LTE populations the equation (|2.2p holds and that for the case of 
LTE, &,- = 1. 



4.2 Radiative rates 



The e xpressions for radiative rates shown here closely follow the book of lMihalad 
(ll97J^ . 



4.2.1 Photoionization and photorecombination transitions 

Photoionization: The radiative rate for a photoionization transition is obtained 
by expressing the amount of absorbed energy E during photoionization from the 
bound state i to the free state k, E = AirJ^aikiv) dv. The number of photoion- 
izations is obtained by dividing the energy E by hv and integrating from the 
ionization frequency I'o to oo. 






(4.5) 



Photorecombination: Photorecombination is a collisional process, since the 
ion catches the free electron during this process. If the medium is in a thermody- 
namic equilibrium, then from the detailed balance it follows that J^, = B^. The 
number of recombinations has its equilibrium value 



(riiRik)* = {rikRki)* = "*47r 



liv 



B„ dv. 



(4.6) 



By adding and subtracting the expression B,ye ''T under the integral and using 
the expression for the Planck function p.4[) , we obtain 



(nkRki)* = <47r 



liv 



2hv^ 



B, 



' I'T Ah 



(4.7) 



Dividing (|4.7p by n^ yields an expression for a rate per ion (the probability). This 
expression has to be valid also outside LTE, so we replace the Planck function B^ 
by actual mean radiation intensity J^- Finally, multiplying the result by actual 
population of the free state Uk we obtain 



rikRki = nk 



An 



aik{v) 
hv 



2hu^ 



+ J^ 



e ''T Av 



(4.8) 



where n*/n*f. — ne$ifc(T) is the so called Saha-Boltzmann factor. 
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4.2.2 Bound-bound transitions 

We consider transitions from the level i to the level /. The number of transitions 
caused by intensity J{v) in the interval dvdLo is niBii(j)ii(i')J^ Av, where (puiv) is 
the line profile function for the given transition, and Bu is the Einstein coefficient 
for absorption. Since aii[i') = a°;(^(i/) = j^Biicf)ii{v)^ we may express the total 
number of absorptions obtained by integration over the line profile as 

riiRii ^TiiBii / 4>ii{v)JijAv ^niA-K / -^ J^dv. (4.9) 

The total number of emissions is given by a sum of spontaneous and stimulated 
emissions 

hv 

(4.10) 

Using the Boltzmann equation (j2.2p . we may write the total number of de-excita- 
tions as 

niRu^m-^^T: ^^i-^ + jAe-T^dv. (4.11) 

4.2.3 Net radiative bracket 

The net rate for the transition i f-> I may be obtained by subtracting upward and 
downward radiative rates 



niRii = ni ( An + Bn I (t)ri{v)J^ dv 1 = ?i;- ^a°, 

hvii gi 



-J,, d, 



ni lAii + Bii / (t)ii{v)Jydv\ -UiBii j (j)ii{v)J^dv = njAuZu. (4.12) 

We have introduced the net radiative bracket Zu as 

Zu = 1 -T-^ / 4>u[iy)Ji, di^ = 1 - ^ , (4.13) 

where Su is the source function (independent of ly) of the transition i <r^ I. For 
the case of detailed radiative balance Zu = 0. 

4.3 Collision a I rates 

Inelastic collisions transfer energy between the thermal particle motion energy and 
the internal energy of atoms. Although collisions between all particles occur in 
stellar atmospheres, it is mostly sufficient to consider only collisions with electrons, 
because the ratio between electron and ion thermal velocities is Wth,c/i'th,i ~ 43'\AA, 
consequently, the ion contribution to collisions can be neglected. 

The upward coUisional rate from the level i to the level I is obtained by inte- 
gration of the total cross section of the transition aii{v) over the electron velocity 
distribution f{v), which is Maxwellian in our case, 

/"OO 

Uidi ^ UiUe crii{v)f{v)vdv = niUequiT), (4-14) 
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wher e the quantity qii{T) is called the effective collision strength (see, e.g. lMihalas . 
Il978l Eq. 5-68). The downward rate follows from the detailed balance in thermo- 
dynamic equilibrium, n'^Cu = n*Cii, and is 

niCu = ni (^) Ca = m (^\ neqa(T). (4.15) 

Thus, for calculation of both upward and downward collisional rates we need to 
know only the effective collision strength qu [T] . 

Detail s of calculations of collisional rat es can be found elsewhere in these pro- 
ceedings (|Butleii . I2OIOI: iMashonkinal . |2010|) . 



4.4 Set of statistical equilibrium equations 

Let us consider an element having NL levels. For each level i we may write an 
equation 

NL NL 

rii J2 iR^l + C^l) + J2^l (Ri- + ^'0 = 0, (4.16) 

1=1 l^i 

l^i 1=1 

where Ra and C^i are given by equations (|43)) . (|4?8l) . (|4J)) . (|4.11l) . (|4.14l) . and 
(|4.15p . Since the system of these equations for all levels is linearly dependent (for 
NL levels only NL — 1 levels are linearly independent), we must supplement the 
equations by some additional equation. 

For model atmosphere calculations, we have two equivalent options (both phys- 
ically and numerically). We can use either the charge conservation equation. 



where Njk means the number density of ions j of the element fc, or the particle 
conservation, 

^5]7V,fe = (iV-ne), (4.18) 

k j 

where A^ is the total number density in the atmosphere, are being used. 

For the case of the solution of the NLTE problem for the trace element, the 
abundance equation 

E^^'^^— E^^^-f (4.19) 

ttref ^-^ 
J 3 

can be used. Here ak is the abundance of the element k and ttref is the abundance 
of the reference atom, which is usually (but not necessarily) hydrogen. Equation 
(|4.19l) relates the abundances of the studied and reference elements. 
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4.5 Non-explicit levels 

Usually the atoms that exist in reality have much more complicated structure than 
we are able to include in our models. As a consequence, some energy levels are not 
explicitly considered in the statistical equilibrium equations. This mainly concerns 
highly excited levels. However, it does not mean that they do not exist and we have 
to include them in the system of the statistical equilibrium equations somehow. 
There are several possibilities of how to handle this task. We may, of course, 
neglect them. However, it is better to include them in the total number density 
and to assume that they are in LTE with respect to the next higher ion. Also the 
corresponding collisional excitation rates to these levels can then be included in 
the generalised (or modif ied) collisional ionization rate (for a detailed description 



of this approximation see lHubenvl . Il988l : iHubenv et al.l . 119941 ). The radiative rates 



to these levels are assumed to be in detailed radiative balance and they do not 
enter into the equations of statistical equilibrium. 

Alternatively, they may be merged into one or more superlevels (see lHubenv fc Lana 

19951) , which can then be included in the set of energy levels and in the equations 



of statistical equilibrium to allow their departures from the local thermodynamic 
equilibrium. 

5 Conclusions 

In this review we summarized the basic equations for the solution of the NLTE line 
formation problem with a focus on the solution of the problem for trace elements in 
stellar atmospheres. It has to be emphasized that in NLTE calculations for trace 
elements preference should be given to background NLTE model atmospheres, 
since using LTE model atmospheres is inconsistent with the condition of NLTE. 
It must be carefully decided which elements can be treated as a trace element and 
for which such treatment is questionable. After obtaining the numerical results all 
assumptions have to be verified, especially that of the traceness. 

The author thanks Drs. Barry Smalley and Glenn Wahlgren for their comments to the manuscript. 
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